Thermoelectric and Magnetothermoelectric Transport Measurements of Graphene 
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The conductance and thermoelectric power (TEP) of graphene is simultaneously measured using 
microfabricated heater and thermometer electrodes. The sign of the TEP changes across the charge 
neutrality point as the majority carrier density switches from electron to hole. The gate dependent 
conductance and TEP exhibit a quantitative agreement with the semiclassical Mott relation. In 
the quantum Hall regime at high magnetic field, quantized thermopower and Nernst signals are 
observed and are also in agreement with the generalized Mott relation, except for strong deviations 
near the charge neutrality point. 



Recent research on graphene, a zero-mass, zero-gap, 
ambipolar material system, has focused on its unique 
electronic properties Although there have been theo- 
retical investigations of thermal and thermoelectric prop- 
erties of this 2-dimensional material [l, 0, 0, [sl, 0, 0, Q , 
only an indirect measurement of the thermal conduc- 
tivity employing Raman spectroscopy has been re- 
ported recently, and thermoelectric transport has yet to 
be explored experimentally. In general, the thermoelec- 
tric power (TEP) is of great interest in understanding 
electronic transport due to its extreme sensitivity to the 
particle- hole asymmetry of a system [lol, [HI • Measure- 
ment of the thermoelectric properties of graphene can 
thus elucidate details of the electronic structure of the 
ambipolar nature of graphene that cannot be probed by 
conductance measurements alone. Furthermore, employ- 
ing the semiclassical Mott relation [10], one can ver- 
ify the validity of the Boltzmann approach to trans- 
port in graphene by comparing the simultaneously mea- 
sured conductance and TEP at different chemical poten- 
tials [ll| . The application of a magnetic field in addition 
to a thermal gradient provides a valuable experimental 
tool to investigate magnetothermoelectric effects, such 
as the Nernst effect in 2-dimensional (2D) electronic sys- 
tems in the quantum Hall (QH) regime |l3l. Il4j. ll5L llq. 

In this letter, we first discuss the temperature and car- 
rier density dependence of the TEP measured in meso- 
scopic graphene samples, and its relation with the Mott 
formula. We then discuss the magnetothermopower mea- 
surement where we observe unusual quantized TEP in the 
presence of the QH effect. 

Graphene samples are fabricated using the mechanical 
exfoliation method on 300 nm Si02/Si substrate. The 
graphene carrier density is adjusted with applied gate 
voltage to the degenerately doped Si substrate. Elec- 
trodes are defined with standard electron beam lithog- 
raphy, followed by Ti/Au (3/30 nm) evaporation and 
liftoff processes. The TEP measurement technique was 
described in detail elsewhere pjj]. In brief, a controlled 
temperature gradient, AT, is applied to the sample by 
a microfabricated heater while the resulting thermovolt- 
age, A V, is measured by the voltage probes to acquire the 




FIG. 1: (a) Conductivity and (b) TEP of a graphene sample 
as a function of Vg for T = 300 K (square), 150 K (circle), 
80 K (up triangle), 40 K (down triangle), and 10 K (diamond). 
(Upper Inset) SEM image of a typical device, the scale bar is 
2/im. (Lower Inset) TEP values taken at Vg = —30 V (circle) 
and -5 V (square). Dashed lines are linear fits to the data. 
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(see the upper inset of Fig. 1 for device 



layout). Local temperature variations are measured with 
two metal 4-probe micro-thermometers. Oxygen plasma 
etching is used to isolate the graphene flake electrically 
from the heater electrode. We use Raman spectroscopy 
and QH measurements in order to verify that our samples 
are single layer graphene [l3, [3 • Devices are annealed 
in a vacuum cryostat to 400 K in order to improve mobil- 
ity and electron-hole symmetry. We measure both the 2- 
terminal mesoscopic conductance, G, utilizing the outer 
electrodes and the 4-terminal conductivity, a, employing 
the Hall bar geometry. To measure the TEP we apply a 
low frequency AC signal to the heater electrode and pick 
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up the resulting 2cj thermovoltage reading from the outer 
electrodes as we change the graphene carrier density with 
the back gate. To stay in the linear response regime the 
condition AT <C T is always satisfied. 

In this work, we study 10 different graphene sam- 
ples with mobility ranging from 1,000 to 7,000 cm^/Vs. 
Fig. la-b show simultaneously measured a and TEP as 
a function of applied gate voltage Vg for a representative 
graphene device in the temperature range of 10 - 300 K. 
As it was shown previously |l|, 0], a oc Vg away from 
the charge neutrality point (CNP). The conductivity be- 
comes minimum at the CNP, corresponding to = Vd, 
where Vd ^ V for this particular device. The peak 
value of TEP reaches ~ 80/iV/K at room temperature in 
the majority of devices measured. Notably, the sign of 
the TEP, which indicates the sign of the majority charge 
carrier, changes from positive to negative as Vg crosses 
the CNP. The electron-hole asymmetry observed in the 
conductance and also in the TEP is most likely due to 
charged impurities in the oxide substrate 1^, 20] or the 
formation of a p-n junction at the electrode contacts [ilj . 
Universal conductance fluctuations (UCF) are observed 
in the conductance for T < 30 K and are quantitatively 
correlated to the oscillations in the TEP, as we will show 
later. 

While a decreases as T increases, the magnitude of 
the TEP increases monotonically as T increases. The in- 
set of Fig. lb shows a linear temperature dependence of 



n (lO^^cm'^) 



TEP at two representative gate voltages, Vg 



-30 V 



and Vg = —5 V. The linearity of the TEP in T sug- 
gests that the mechanism for thermoelectric generation 
is diffusive thermopower. Phono n-drag components of 
the thermopower, which usually appear in the form of a 
non-monotonic temperature dependence of TEP [22], are 
not present here, which is consistent with weak electron- 
phonon coupling in graphene 2^, 13]. 

Within the Boltzmann formulation [25], the Seebeck 
coefficient, 5, equivalent to the TEP, represents the en- 
tropy transported per unit charge. This thermodynamic 
relation provides us with a qualitative understanding of 
the overall behavior of the observed S{Vg) as Vg tunes 
the carrier density, n, and thus the Fermi energy, Ef, 
of the graphene samples. Far away from the CNP, the 
electronic system is highly degenerate and the entropy 
transported in the channel is proportional to the number 
of thermally activated carriers over the degenerate Fermi 
sea, leading to 5* ~ T/Ep. Close to the charge neutrality 
point, however, the electron gas becomes non-degenerate 
and both electrons and holes are thermally populated. 
Here the TEP gains opposite contributions from electrons 
and holes, so l^l decreases as |n| decreases. However, the 
non-degenerate limit is difficult to achieve in most of the 
experimental temperature range since the electron-hole 
puddle formation near the charge neutrality point yields 
a finite local density >10^^ cm~^ [^26], corresponding to 
^100 K of energy broadening. 
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FIG. 2: Measured TEP values (thick red) and TEP values 
calculated with the Mott formula (thin black) using Eq. (1) 
as a function of applied gate voltage for T = 200 K (a), 40 K 
(b), and 15 K (c). 



For a quantitative understanding of the TEP data, we 



now employ the semiclassical Mott relation p^, Il2[ . In 
mesoscopic graphene samples, this relation can be tested 
by tuning Ep using Vg'. 
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3|e| 



G dVg dE 



(1) 



E=Ef 



where e is the electron charge and /cb is the Boltzmann 
constant. This formula can be evaluated by numerically 
differentiating the measured G with respect to Vg and 



considering ^{Vg) 



VlM^, where /\Vg 



Vg — Vdj ^ is Planck's constant, the Fermi velocity vp = 
lO^m/s, and the gate capacitaince Cg = 115 aF//im^ 
in our device geometry [18]. A comparison between the 
measured TEP and that predicted by the Mott formula is 
plotted in Fig. 2 for temperatures T =15, 40, and 200 K, 
where we find excellent agreement at lower temperatures. 
We note that Eq. (1) is in the form of the mesoscopic 
formula \2l\ w here the conductance includes the contact 
resistance [28]. At higher temperatures, local thermal 
equilibrium is established, and the mesoscopic form Eq. 
(1) using G becomes less accurate, as it is already seen 
in the 200 K data in Fig. 2. Such deviation grows at 
higher temperatures. Recent theoretical work [7] predicts 
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FIG. 3: (a) Hall conductance; (b) longitudinal TEP Sxx'i 
and (c) transverse TEP Syx as a function of Vg and B at 
T = 10 K. The insets shows the electrode configurations for 
Sxx and Syx measurement, respectively, (d) Rxx; (e) Rxy; (f) 
Sxx'i ctnd (g) Syx as a function of Vg at a fixed magnetic field 
5 = 8.8 T and T = 4.2 K. The vertical doted lines indicate 
filling factors corresponding to zy = ±2, ±6, ±10, ±14. All 
the data is averaged with respect to positive and negative 
values of magnetic field in order to remove mixing between 
longitudinal and transverse components. 



that in graphene the enhanced inelastic scattering time of 
carriers in the presence of electron-electron interactions 
yields TEP that deviates strongly from the Mott relation. 
Further experimental work with higher mobility samples 
(in current samples disorder dominates transport even 
at room temperature) is needed to elucidate this non- 
degenerate, high temperature limit of TEP in graphene. 

We now turn our attention to the magneto-TEP. In 
the presence of a perpendicular magnetic field, a Lorentz 
force bends the trajectories of the thermally diffusing car- 
riers to produce the Nernst TEP, i.e., transverse com- 
ponents of the TEP, Syx' For this measurement, we 
utilize a multiple probe configuration to measure the 
tensor components of resistivity Rxx and Rxy and also 
TEP components Sxx and Syx as shown in the inset 
of Fig. 3(b-c). Fig. 3a displays the Hall conductivity 
o'xy = ~ Rxy /[Rxy + (WRxx/L)'^] as a fuuctiou of mag- 
netic field and gate voltage, where W and L are the width 
and length of the Hall bar device, respectively. The lin- 
ear growth of the oscillations of cFxy{Vg) with increasing 
magnetic field clearly indicates that Shubnikov-de Hass 



oscillations (SdH) are present in graphene as previously 
observed [13] • Remarkably, similar oscillatory fea- 
tures develop in both Sxx and Syx as shown in the Fig. 
3(b-c), suggesting that magneto-TEP also provides a sen- 
sitive probe for the SdH effect in graphene. 

In the high field regime, 5 ~ 9 T, fully developed QH 
effect is observed for filling factors u = ±2, ±6, ±10, ±14, 
indicated by the quantized plateaus in Rxy and zeros in 
Rxx- Vertical dashed lines in Fig. 3(d-g) indicate bound- 
aries between LLs where we find Sxx and Syx tend to 
zero since no carriers are available to participate in dif- 
fusion. At the Landau level centers \Sxx\ is maximal and 
as higher LL are occupied these peaks decrease. Theoret- 
ically, this quantized behavior of Sxx is expected in the 
QH regime of 2D electron gas systems [13|, ll6|. In the 
low temperature limit, the peaks of Sxx are predicted to 
be given by S'!^J^^ ~ ^^"-^^2 clean samples, where 
N in graphene indicates the LL index. Experimentally, 
thermal and disorder induced broadening of LLs tends to 
broaden Sxx^ making the quantized values smaller than 
the above prediction 15|, ll6| . The observed S'^^^ shows 
the quantized 1/N trend discussed above, but with a re- 
duced factor, indicating the presence of disorder in our 
graphene samples. 

Further analysis which includes the Hall conductivity 
tensor and the Nernst TEP is performed by employing 
the generalized Mott formula that holds for all tempera- 
tures and can be applied to the QH regime pj, : 
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where a^j is the conductivity tensor and z,j represent 
the X and y components. According to Eq. (2), we ex- 
pect that Sxx/T becomes temperature independent in 
the QH limit where ciij is temperature independent. In- 
deed, Fig. 4a shows that the measured Sxx/T is rather 
temperature independent for T < 20 K. Above this tem- 
perature, however, the oscillation amplitude of Sxx/T 
becomes smaller since thermal activation across the dis- 
order broadened LLs becomes appreciable and the QH 
effects disappear. Fig. 4b shows good agreement between 
the measured Sxx{Syx) and the calculated S^J'^^S^J'^^) , 
indicating that the semiclassical Mott relation, extended 
into the QH regime, works well for both holes {N < 0) 
and electrons {N > 0). Near the zeroth LL {N = 0) at 
the CNP (Vg ^ Vd), however, the measured Sxx exhibits 
a pair of anomalous oscillation that shows a distinct de- 
viation from Eq. (2). We remark that this pair of Sxx 
peaks near the zeroth Landau level shows opposite TEP 
polarity, i.e., positive peak on the electron side {Vg > Vd) 
and negative peak on the hole side {Vg < Vd)- This 
pair of opposite polarity oscillations could be indicative 
of a peculiar QH state stemming from the = LL. 
It has been speculated that this LL produces a pair of 
counter propagating edge states when the Fermi level is 
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manuscript, we became aware of a related work 
with similar experimental results [30i] . 



FIG. 4: (a) Temperature dependence of the Sxx/T as a func- 
tion of Vg at B = 8.8 T. (b) Measured Sxx (thick red) and the 
Sxx predicted by the generalized Mott relation (thin black) 
at 5 = 8.8 T and T — 4.2 K showing diviations near the 



iV = LL. (Inset) Corresponding S., 



Measured 



and S, 



Mott 



at the CNP [29]. In this scenario, the polarity of the 
inner edge state is opposite to the sign of the bulk ma- 
jority carrier, which may generate the anomalous TEP 
observed. In addition, a large enhancement of ^^ya; (Inset 
of Fig. 4b) is observed near the CNP, which also deviates 
strongly from S^^^^. Such enhancement of the Nernst 
signal is predicted in conventional 2D electron systems, 
where the magnitude of Syx strongly depends on the dis- 
order strength [13, [3] • 

In conclusion, we report carrier density depen- 
dent measurements of the thermoelectric coefficient in 
graphene which scales linearly with temperature. The 
semiclassical Mott relation explains the observed carrier 
density modulated TEP particularly well in the low tem- 
perature regime. In the QH regime at high magnetic 
field, the tensor components of the TEP show quantized 
behavior along with unusual deviations from the gener- 
alized Mott relation near the charge neutrality point. 
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